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We study sustained oscillations in two-dimensional oscillator systems driven by Rayleigh-typc 
negative friction. In particular we investigate the influence of mismatch of the two frequencies. 
Further we study the influence of external noise and nonlincarity of the conservative forces. Our 
consideration is restricted to the case that the driving is rather weak and that the forces show 
only weak deviations from radial symmetry. For this case we provide results for the attractors and 
the bifurcations of the system. We show that for rational relations of the frequencies the system 
develops several rotational excitations with right/left symmetry, corresponding to limit cycles in 
the four-dimensional phase space. The corresponding noisy distributions have the form of hoops 
or tires in the four-dimensional space. For irrational frequency relations, as well as for increasing 
strength of driving or noise the periodic excitations are replaced by chaotic oscillations. 
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I. INTRODUCTION 



The first comprehensive theory of nonlinear oscilla- 
tions was developed by Lord Rayleigh in the y ears 1883- 
1894 and is represented in his pioneering book ijRavleighl 
1894). The present state of art of nonlinear dynamics 
was deeply in fluenced by the pioneering work of Leonid 
Shilnikov (see IShilnikovi Il997t IShilnikov et~aH Il998albl 
and refe rences therein) whi c h is reflected e.g. in th e 
books of lAnishchenkcl (|l995[klAnishchenko et al\ (|2002^ . 
In the last decade many investigations were devoted to 
the st ochastic dynamics of dissipative hamiltonian sys- 
tems ijAnishchenko et adl2003lKlimontovichlll994|l . 

In this paper we extend previous studies on rota- 
tional excitations of two-dimensi onal oscillators driven 
by Ra yleigh-type negative friction ijErdmann et q^ . l2002l 
2000). The study of driven oscillatory modes on a plane 
has interesting applications for modeling animal mobil- 
ity all the way from microorganisms like bacteria and 
Dictyostelium discoedium slime mold llBen- Jacob et all 
l200(ilCzir6k et adll996HRappel et oZl.ll999l) upto flock s 



of birds l|Toner and TuLll993: Weimerskirch et artl200ll) 
schools of f ish ( Hubbard_&LaUl2004; Inad a and Kawachl 
l2002HNiwal 11994 iParrish et adl2002D. swarm s of Da ph- 
nia l)Erdmann et all 2004 lOrdemann et all 12003^ or 
wildebeest ijTopaz anc . Bertozzil |2004|) . For example we 
introduced in earlier work the general idea of stochasti- 
cally moving species, active Brownian particles. We want 
to recall this approach which will be used later on. Ac- 
tive Brownian particles are Brownian particles with the 
ability to take up energy from the environment and use it 
for the acceleration of motion. Simple models composed 
of active Brownian particles were studied in many ear- 



lier works (e.g.lHelbing and Molnailll995tlKlimontovichl 
11994 ISchienbein and Gruleit Il993ft . 



A specific problem we would like to address here 
is: what is the consequence of broken radial symme- 
try (asymmetry of the conservative forces). In par- 
ticular we are interested in the problem: Can sponta- 
neous rotations be stopped by certain a mount of asym- 
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SteuernageleT^I 11994) which are able to convert stored 
internal energy into motion. In this paper, we will study 
only the motion, especially rotational ones, in external 
fields on a two-dimensional plane. 



The article is organized as follows. In Section [n] we 
introduce the equations of motion, the pumping by neg- 
ative friction and outline the basic dynamics of our model 
including Langcvin equations. In Section IIIII we out- 
line our previous studies of rotationally symmetric ex- 
ternal potentials and discuss the attractors for rational 
frequency relations. In Section llVI we give an analysis of 
frequency mismatch and of Arnold-type bifurcations of 
linear oscillators as a function of the mismatch. Further 
we study in Section [3 the limit cycle attractors of non- 
linear oscillators with radial symmetry. In Section IVII 
we investigate the influence of noise on the attractors 
and show that noise leads to a broadening on the line- 
attractors. 
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II. DYNAMIC EQUATIONS FOR TWO-DIMENSIONAL 
OSCILLATORS 



III. ATTRACTORS FOR LINEAR OSCILLATORS DRIVEN 
BY RAYLEIGH FRICTION 



The dynamics of the systems studied here is based on 
Langevin equations, kno wn from the theory of conven- 
tiona l Brownian motion ijHanggi et all . ll99Ct iLangevinl 
1908). For the two-dimensional space we get four first 
order coupled differential equations in the phase space 
{xi,x 2 ,Vi,v 2 }: 
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In this dynamics we assumed three kinds of forces in the 
dynamics: 

1. Conservative external forces generated by the po- 
tentials U (xi, x 2 ), 

2. nonlinear dissipative forces modeled by the friction 
(-«7(«)), 

3. stochastic forces assigned by £(t). 

Several linear and nonlinear conservative forces will be 
introduced and studied in the following sections. The 
dissipative forces are modeled by the friction function 
~f(v). This function is the source of dissipative interac- 
tions with the surrounding. In equilibrium the friction is 
passive 



7(f) = 7o = const. 



(2) 



We will consider here in more deta il active f r iction mod- 
eled by the classical Rayleigh law ijRavleighl Il894|) : 



70) = 7(^1^2) =a-0(vf + v 2 ). 



(3) 



The stochastic forces are modeled by white Gaussian 
noise £ with vanishing mean and 



mw)) = 5(t-t') 



(4) 



and scaled with strength D . In equi l ibrium and in case 
of passive friction following lEinsteinl l)l956fl one gets an 
energy balance between the strength of the stochastic 
force, D, and the passive friction acting on the object. It 
is expressed by the simple fluctuation-dissipation relation 



D 



70 



(5) 



where 9 = ksT is a measure for the temperature. 



We will study in this section two-dimensional linear 
oscillators described by the potential 



m 



U{xi,X 2 ) = — (WjZ! +u 2 x 2 ) 



(6) 



which are driven by Rayleigh-typ e negative fricti on as 
in Eq. (j3J| . It is well known since iRavleighl l)l894j) that 
in the one-dimensional case the system possesses a limit 
cycle corresponding to sustained oscillations with the en- 
ergy Eq = a/ 8. The two-dim e nsion al case is much more 
complicated. lErdmann et all l)2000|) have shown for the 
symmetrical case with oj 2 = lo\, that two limit cycles 
in the four-dimensional phase space are developed. The 
projection of these sustained oscillations on the {vi,v 2 }- 
plane and on the {xi, x 2 }-p\a,ne are circles 



2 _,_ 2 

X 1 -|- X<2 



o 



The limit cycle energy is 



En = 



(7a) 
(7b) 

(8) 



lEbeling et all |)1999|) have shown, that any initial value 
of the energy converges (at least in the limit of strong 
pumping) to 



H 



Eq = mw, 







(9) 



This corresponds to an equal distribution between kinetic 
and potential energy i.e. both parts contribute the same 
amount to the total energy. The motion on the limit 
cycle in the four-dimensional space may be represented 
by the four equations 

xi = tq cos(woi + 4>o) (10a) 

vi = — tqloq sin(w * + 4>o) (10b) 

x 2 = r sin(w t + 0o) (10 c ) 

v 2 = r uJo cos(uj t + 0o ) (lOd) 

The angular frequency follows by estimations of the time 
the particle needs for one period moving on the circle of 
radius ro with constant speed vq: 



vo 

w = — 
ro 



UJl = UJ 2 



(11) 



This means, the particle rotates even at strong pump- 
ing with the frequency given by the linear oscillator fre- 
quency uji — u) 2 . The trajectory defined by Eqs. l|TU|) 
is an exact solution of the dynamic equations describ- 
ing the first (forward) limit cycle. The shape of the 
trajectory is like a hoop in the four-dimensional space. 
Most projections to the two-dimensional subspaces are 
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circles or ellipses however there are two subspaces namely 
xi,z>2} and {x^fil w here the projection is like a rod 
Erdmann et all l200fl) . Reversing the initial velocities 
of the system, a second limit cycle can be obtained. This 
limit cycle forms also a hula hoop which is different from 
the first one. However both limit cycles have the same 
projections on the {xx,X2} and on the {vi, i>2}-plane. 
The projection to the {aq, a^j-plane has the opposite 
sense of rotation in comparison with the first limit cy- 
cle. The projections of the two hoops on the {x\, X2}- 
plane or on the {v%, i^j-pl&ne are two-dimensional rings. 
The hoops intersect perpendicular {xi, V2}- and {X2, v±}- 
planes. The projections to these planes are rod-like 
and the intersection manifold with these planes consists 
of two ellipses located in the diagonals of the planes 
l|Erdmann et aZJ.l2000^ . 

So far we repeated known results for the symmetrical 
case. Dynamical systems with radial symmetry are de- 
generate and structurally unstable in the mathematical 
context. From the physical viewpoint, radial symmetry 
is a special situation, i.e. the gravitational field of point 
masses or the Coulomb field for charges has strict radial 
symmetry. Therefore radial symmetry holds also for a 
two-dimensional mass-point pendulum. In real physical 
systems the radial symmetry is in general broken, e.g. 
a real pendulum in the earth field has no strict radial 
symmetry. Thus the oscillator with radial symmetry can 
only be considered as a particular case of correspond- 
ing real system which often has some asymmetry. Thus 
our motivation is, to investigate the consequences of bro- 
ken radial symmetry and deviations from linearity on the 
generation of oscillatory modes. 

First we are going to study systems with frequency 



mismatch u>2 7^ u)\. In the conservative case 7 = and 
D = we find a dense set of exact solution for the deter- 
ministic dynamics given by 

Xl (t) = Ax cos(wii + $1) (12a) 

Vl (t) = -Axujxsin(ujxt + $x), (12b) 

af 2 (*) = A 2 sxxi(Lu 2 t + <f>2) (12c) 

v 2 {t) = A 2 uj 2 cos(Lu 2 t + $2) , (12d) 



where A4 and $j are given by the (arbitrary) initial con- 
ditions. Let us study now the driven case with negative 
friction according to the Rayleigh law without noise. 
Then for a > the above periodic solution 1)1 2H would 
remain to be a solution, if 

2,2 2 a r-t o-v 

vx +v 2 =v Q = - (13) 

is fulfilled. In the case of symmetrical oscillators w\ = 
u>2 = u>o we could fulfill this condition in an exact way 
by a special choice of the amplitudes: 

Ax = A 2 = r = ^ (14) 

LUq 

This means, we find instead of a dense set of exact so- 
lutions just one attracting solution (which still is ex- 
act). This periodic and stable exact solution represents 



a limit cycle corresponding to a circular path even at 
strong pumping with the frequency given by the simple 
harmonic oscillator frequency loq. However in the case 
u>x ll>2 the situation is much more difficult than in the 
symmetrical case. First we are looking for approximate 
stable solutions in the case of rational relations of the fre- 
quencies 0J2 = nujQ and u>x = rriLUo with n,m — 1,2,.... 
The condition <|13|l is in average over one period fulfilled, 
if 

A\u\ = A\u\ = v 2 Q (15) 
This way we get the stable amplitudes and phases 

Ax = — ; A 2 = — (16a) 
lox ^2 

*i = -|; $2 = +f ( 16b ) 

We will show that these approximative solutions de- 
scribe at least for small relations n : m again a pair of 
forward/backward limit cycles. By introducing the ap- 
proximation ^| into Eq. i|12|) we find for the n > m the 
following analytical approximation for the limit cycles: 



Xx{t) = ro sin(ma;oi) (l^ a ) 

Vx(t) = romojQ cos(mu>ot) (17b) 

X2(t) = rocos(nu>ot) (17c) 

V2(t) — — ronojQ $m.(nu>Qt) . (17d) 



The curves obtained by this approximation for n : 
m = 2,3,4,5 are shown in Fig. ^ I n the analyti- 
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FIG. 1 Projections of asymmetric limit cycles to the {xi, X2}- 
plane for several cases of rational relations between the two 
frequencies m : n = 2, 3, 4, 5 (analytical approximation for the 
attractors (Eqs. il'Il and 1161 1. 

cal approximation given above, the forward limit cycles 
and the backward limit cycles have identical projections 
on the {xx, a^j-plane. However our analytical formula 
(Eqs. (|12|) and (|16|l 1 is only a rough approximation, as 
we will demonstrate by simulations. We show in Figs. [21 
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FIG. 2 Projections of the two limit cycles to the {x±,X2}- 
plane corresponding to m : n = 2 resonance obtained from 
simulations (Rayleigh law: a = 5, /3 = 1). 
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FIG. 3 Projections of the two limit cycles to the {«i,«2}- 
plane corresponding to m : n — 2 to the {«i , V2 }-plane ob- 
tained from simulations (same parameters as in Fig. |5J . 



and |3 the result of simulations for the two attractors in 
the case m : n = 2. Clearly, the two attractors have 
different projections on the {x\, ar2}-plane (see Fig. 0) 
as well as on the {vi, «2}-plane (see Fig. EJ. As we see, 
the projection differ in amplitude and phase a little bit 
from the analytical approximation; in particular clock- 
wise and counterclockwise limit cycles are shifted in the 
projections. It is interesting to note that the projections 
to the {xi, v i}-plane and to the {x2, «2}-plane shown in 
Fig. 0| are equal for the clockwise and the counterclock- 
wise limit cycles. 

The rather complex winding structure of the attractor 
for n = 2 in the four-dimensional phase space can be 
guessed by looking at projections on three-dimensional 
subspaces as demonstrated in Fig. [3] 

The results obtained from simulations for n = 3 are 
shown in Figs. El and We see that the attractors ob- 
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x 

FIG. 4 Projections of the two limit cycles to the {0:2, V2}- 
plane for m : n = 2 and respectively obtained from simula- 
tions (same parameters as in Fig. |5J . Clockwise and counter- 
clockwise limit cycle have practically the same projections. 




FIG. 5 Projections of the two limit cycles to the {xi, X2, V2}- 
plane corresponding to the m : n — 2- resonance obtained from 
simulations (same parameters as in Fig.|5J. 



tained from simulations are more complex than the from 
the nice and rather symmetrical curves obtained from 
the analytical approximation for the case n = 3 which 
were presented in Fig. ^ With increasing rational n the 
attractor fills more or less dense a rectangular region. 
For irrational values of W2/W1 the trajectories are always 
dense in a nearly rectangular region of the coordinate 
space. This case will be studied below. 



IV. BIFURCATION ANALYSIS OF THE LINEAR 
OSCILLATORS WITH FREQUENCY MISMATCH 

Our study was limited so far to small rational fre- 
quency relations m : n and small or moderate strength of 
driving (small positive values of the bifurcation param- 
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FIG. 6 Projections of the two limit cycles to the {xi,X2}- 
plane corresponding to the m : n = 3-resonance obtained 
from simulations (u>i — 2.7, u>2 = 1, all other parameters as 
inFig.HJ. 




FIG. 7 Projections of the two limit cycles for m : n — 3 to 
the {xi, X2, «2}-plane obtained from simulations. 



eter a). In order to introduce the consequences of a an 
irrational frequency mismatch between the partial sub- 
systems in a more general setting we write the expression 
for the potential Eq. (JSJ as follows: 



U(xi,x 2 ) 



2/ 2 , \2 2\ 
■UJ (X 1 + A X 2 ), 



A = ^ 

0J 2 



CJ = OJl- 



(18a) 
(18b) 



part of Eq. JQ may be written as 

Xi = Vi (19a) 

ill = [a - (vf + v$)] vi - u$xi (19b) 

±2 = v 2 (19c) 

v-2 = [a-p{v{ «a -^o A ^2 (19d) 

The system (|19|l is structurally stable or is one of the 
commo n propos i tions according to V. I. Arnold's nomen- 
clature llArnoldllT96l . It can be interpreted as a model 
for one oscillator on the plane or as a model for two in- 
teracting linear oscillators. 

To understand the dynamics of system (|19|l in the gen- 
eral case, we introduce a complex variable z = x\ + jx 2 . 
Let us first turn again to the symmetric case where 
u>i = uj 2 = u>q- Then from 1|19|) follows 



(20) 



z-/?(^--|i| 2 Jz + ^ 2 z = 0. 
Equation (|20|) has periodic solutions of the form: 



z(t) = zexp(±jw t) = \z\ exp(j$) exp(±juj t), (21) 

where the phase $ takes any value in the interval [0, 2ir]. 
When we consider the symmetric case, we have an infinite 
number of periodic solutions (see Eq. However, 
linear analysis cannot yield information about their sta- 
bility. In numeric calculations, we can detect several limit 
cycles each of them po ssessing its own type of symmetry. 
As bv lErdmann et al. (2002) has been investigated, the 
stability of periodic solutions can can be predicted by the 
calculation of the Floquet multipliers of the fixed points 
within the Poincare section of the limit cycle. When there 
is a detuning (A ^ 1), only two limit cycles, representing 
clockwise and counterclockw ise rotations remain s table 
as has already been found bv lErdmann et all l|2000|h All 
other solutions have at least one multiplier leaving the 
unit circle on the complex plane. The two limit cycles 
representing the rotations of a particle in the coordinate 
space are situated within the first Arnold tongue whereas 
outside this region the existence of the cycles vanishes 
and the particles move irregularly. Increasing the de- 
tuning A any further new stability regions appear which 
represent the period-n cycles of the oscillators as have 
been shown in Figs. 1 1171 The parameter region where 
stability of periodic motion can be observed is shown in 
Fig. El 



V. NONLINEAR OSCILLATORS WITH RADIAL 
SYMMETRY 



In the general case the potential U(xi,x 2 ) has an elliptic, 
slightly extended shape, the asymmetry is measured by 
the parameter A. For irrational values of A the solutions 
in the conservative case a = /3 = fill densely a box in the 
phase space. In the general driven case the deterministic 



In the present section we will discuss several exten- 
sions of the theory developed in the previous section to 
nonlinear oscillators. However, for simplicity, we will re- 
strict our study to systems with rotational symmetry, 
U = U(r) with U(0) — and U(f)— monotonically 
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FIG. 8 Stability regions (Arnold tongues) of the m : n res- 
onances. Within the light gray areas the limit cycles as de- 
scribed above are stable. The arrows show the parameter set 
of Figs. HE 



increasing. For the general case of radially symmet- 
ric but anharmonic potentials U (r) the equal distribu- 
tion between potential and kinetic energy is violated. In 



other words the relation 



which leads to 



u>o = Vo/to = w is no more valid. This relation has to be 
replaced by the more general condition that on the limit 
cycle the attracting radial forces are in equilibrium with 
the centrifugal forces. This leads to 



several solutions, the dynamics might be much more 
complicated. An interesting application of the theoret- 
ical results given above is the case of Coulomb forces 
ijSchimanskv-Geier et aZ.Ll2005|) . Another possible appli- 
cation is the following: Let us imagine a system of Brow- 
nian particles which are pairwise bound by a Lennard- 
Jones-like potential U(\r\ — r^\) to dumb-bell-like con- 
figurations. Then the motion consists of two indepen- 
dent parts: The free motion of the center of mass, and 
the relative motion under the influence of the potential. 
As a consequence, the center of mass of the dumb-bell 
will make a driven Brownian motion but in addition the 
dumb-bells are driven to rotate around there center of 
mass. What we observe then is a system of pumped 
Brownian molecules which show driven translations with 
respect to their center of mass. On the other side the in- 
ternal degrees of freedom are also excite d and we observe 
driven rotations. lErdmann et all l)2002|) have shown that 
this approach can be extended to systems of many par- 
ticles in the sense that as far as the mean field of the 
interaction potential can be approximated like potentials 
of the shape of Eq. rotations are going to be sta- 
ble within the first Arnold tongue (Fig. 0. In this way 
we have shown that the mechanisms described here may 
be used also to excite the internal degrees of freedom of 
Brownian molecules. 



VI. THE INFLUENCE OF NOISE 



''0 



If vq is given, the equilibrium radius may be found from 
the implicit relation 



«o 2 = ?V(ro)| 



(23) 



Then the frequency of the limit cycle oscillations is given 
by 



\U'(r )\ 



In the case of linear oscillators this leads us back to pre- 
vious result given in Sec. IIIII For the case of quartic 
oscillators 



Mr) = -r 4 
4 

we get the limit cycle frequency 



wo 



1/2 



(22) The main effect of noise is the spreading of the deter- 

ministic attractors. Let us consider here only the case 
of radially symmetric potentials U(r) which have a mini- 
mum U = at r = and are monotonically increasing 
with r. Then as shown above, the system has two limit 
cycles in the four-dimensional space which are hoop-like 
and have projections at the {xi, X2}-space which are cir- 
cles with the radius tq . This radius is determined by the 
equilibrium between centripetal and centrifugal forces for 
right/left rotations on the circle with the radius tq. In- 
cluding stochastic effects the two hoop-like limit cycles 

(24) 

are converted into distributions looking like two embrac- 
ing hoops with finite size, which for strong noise convert 
into two embracing tires in the four-dimensional phase 
space. 

The transformation of the limit cycles into tires can 
easily be obtained from simulations including white ad- 

(25) ditive noise. Several results for both symmetric (A = 1) 
and asymmetric (A ^ 1) parabolic potentials are demon- 
strated in Fig. [5] The probability distributions may 
be obtained as solutions of the Fokker-P lanck equation 
for th e probability distribution P(r, v, t) ( Kli montovichl 

( 26 ) 1994) 



The explicite solution IjlOjl remains valid, i.e. we find 
again an exact analytical description of the pair of limit 
cycles. For monotonically increasing potentials there ex- 
ist just one stable radius tq. If the equation (|23|) has 



dP 
~dt 



dP 
Or 



VU{r) 



dP 
dv 



(27) 



d ( dP 
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FIG. 9 Distributions according to simulations of 2500 Active Brownian particles on limit cycles in symmetric (a-b) and 
asymmetric (c-f) parabolic potentials. We show projections of simulations on the {xi, £2}-plane and on the subspace {xi, X2, V2}- 
Parameters: a = 5,/3 = 1, D = 0.01. 
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To find explicite solutions of Eg ll2*?l) is a very diffi - 
cult task ijAnishchenko et g;J . l2^3 lKlrnionto vichlll994j) . 
Only in the force- free case 17 = 0, the solution is ele- 
mentary. The stationary solu t ion re ads for the Rayleigh- 



mentary. inc stationary solution rce 
model |Erdmann efaZl . 120021 12000): 



Pq(v) = C exp 



av 
~2D 



(28) 



The shape o f this velocity d i stribu tion Eq. I|28|) can be 
seen in (e.g. lErdmann et all 12002^1 . The bifurcation to 
a limit cycle at the transition from negative to positive 
a occurs for the noisy system as a qualitative change 
of the shape of the distribution from a Maxwell-like to 
a hat-like shape. It is obvious that the system above 
the bifurcation point is far from equilibrium and shows 
a permanent active motion of the particles. Applying 
similar arguments to the stochastic motion in confining 
potentials we expect that the two hoops are converted 
into a distribution with the appearance of two embracing 
hoops with finite size, which for strong noise converts into 
two embracing hoops in th e four-dimensional phase space 
fsee lErdmann et "aZLl2000l for details). In order to obtain 
the explicit form of the distribution, we may introduce 
the amplitude-phase representation 



Vi 
X2 
V2 



pcos(uiot + 4>) 
—pujQ sin(wo£ 
psin(o;ot + <f>) 
pojQ cos ( u>ot + 



(29a) 
(29b) 
(29c) 
(29d) 



crater is determined by the two deterministic limit cycles 
(see Fig. rrU|) . The velocity distribution given by Eq. l(2^|) 
remains to be exact for all radially symmetric potentials. 
The full stationary probability in the four-dimensional 
space has the form of two hula hoop distr i bution s 
ijDeng and Zhul. 12004 ISchimanskv-Geier et all l2005|) . 
The projections of the distribution onto the {£1,2:2}- 
plane and to the {vi, i^j-plane are two-dimensional rings 
(see also Fig. E{a-b)). As in the deterministic case 
the hula hoop distribution intersects perpendicularly the 
{xi, «2}-plane and the {X2, vi}-plane. Again, the pro- 
jections to these planes are rod-like, and the intersec- 
tion manifold with these planes consist s of two ellipses 
located in the diagonals of the plan es ijErdmann et all 
2000: ISchimanskv-Geier et all 120051) . In the determinis- 
tic case one of two the rotational motions within the con- 
fining potential is excited, this rotation remains a stable 
solution of the trajectory of one particle. To this rota- 
tion belongs a certain value of the angular momentum. 
For non-vanishing stochastic perturbations, the particle 
is able to cross the separatrix between the two rotational 
modes (limit cycles). Due to this ability of the particles 
one can observe, sometime s, an inversion of the a ngular 
momentum of the particle ijErdmann et all [2002) . Note 
that one initial condition will be sufficient to rotate cither 



clockwise or counterclockwise (see Fig. 9(b) 



VII. CONCLUSIONS 



where radius p and phase 4> are slow and fast stochastic 
variables respectively. By using the standard procedure 
of averaging with respect to the fast phases we obtain for 
the Raylcigh- model of pumping the following distribution 
of the radii: 



P (p) ~ exp 



D 



2rl 



(30) 




FIG. 10 Probability distribution Po(p) for a radially symmet- 
ric potential in the coordinate space {xi, £2}- The parameters 
a,ro,vo, D are set equal to 1. 

This distribution is (in the present approximation) uni- 
versal and valid for any radially symmetric potential of 
the type specified above. We see that the probability 



We study here two-dimensional oscillations with non- 
linear Rayleigh-type active friction, extending earlier 
work to 



1. oscillations with frequency mismatch 



2. oscillations with radially symmetric nonlinear at- 
tractive forces. 

We investigate several right /left symmetric pairs of limit 
cycles corresponding to 1:1, 2:1 and 3:1 resonances and 
their location in the stability regions (Arnold tongues). 
We show that in the presence of noise the limit cycles 
are converted into hoops in the four-dimensional space 
and give analytical estimates for the probability distri- 
butions of the coordinates and velocities. Thinking of 
applications especially for coherent behavior in animal 
motion on could think of looking for interacting particles 
where the interaction could be approximated (in mean 
field approxima tion) as a devia t ed pa rabolic potential. 
In extention to ijErdmann et all l2002j) on should inves- 
tigate more the question if there would be any coherent 
motion which could be located within the higher Arnold 
tongues. How would the motion of a swarm of particles 
look than? 
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